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Abstract 


-  Homogeneous  and  internally  structured  material  coatings  are  studied  numerically  as 
candidates  for  drag  and  turbulence  generated  noise  reduction  on  submersible  hulls. 
Dynamic  motions  of  compliant  surfaces  may  conceptually  interrupt  random  turbulent 


motions,  reducing  the  turbulent  interface  stresses  and  absorbing  turbulent  energy,  and 
consequently  reducing  the  turbulence  at  the  flow/coating  interface.  Time  dependent,  two 
and  three  dimensional  Monte-Carlo  turbulent  pressure  field  models,  and  interactive 


Navier-Stokes  pseudo-spectral  methods  are  used  to  represent  the  unsteady  flow,  while 
finite  element  methods  are  applied  to  represent  a  variety  of  homogeneous,  layered,  and 
internally  structured  material  coatings.  The  influence  of  added  mass  and  in-depth 
overburden  on  the  material  response  in  ocean  water  is  discussed.  Promising  compliant 
material  coatings  include  sandwiches  of  soft,  homogeneous  layers  between  thin,  stiffer 
elastic  materials  and  internally  structured  coatings  combining  streamwise  ribs,  and 
spanwise  voids  separated  by  stiff  elastic  supporters  imbedded  in  soft  yielding  viscoelastic 


layers. 
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I.  Introduction 


After  the  Second  World  War  interest  in  the  possibilities  of  reducing  drag  through  the 


use  of  compliant  coatings  was  stimulated  by  study  of  the  summary  efficiency  of 


mammals,  particularly  porpoises.  Simulated  skin  coatings  were  developed  and  tested  in 
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the  United  States  by  M.  O.  Kramer  in  the  1950’s  and  1960’s.  ’  More  recently. 


experimental  USSR  results  indicate  that  significant  reductions  in  turbulent  skin  friction 
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drag  are  obtained  through  use  of  thin,  elastically  deforming  surface  coatings. 


Unfortunately,  the  promising  experimental  results,  have  not  been  independently  and 


substantially  reproduced.  The  issue  remains  open. 


Our  present  investigations  focus  on  compliant  coatings  that  exhibit  favorable 


dynamic  response  to  both  laminar  and  turbulent  boundary  la^sr  motions.  Favorable 
response  is  that  associated  with  reduction  of  skin  friction  drag,  in  our  discussion. 


Identification  of  the  fundamental  fluid/solid  dynamic  coupling  mechanisms  having  a 


potential  for  reducing  skin  friction  drag  is  a  primary  goal.  We  emphasize  materials  and 


coatings  whose  characteristics  can  be  altered  by  passive  changes  to  internal  structure,  or 


composite  material  stiffness,  or  dynamic  variations  in  internal  stress,  or  by  wall  heating. 


Such  coatings  should  be  effective  over  a  useful  range  of  flow  conditions,  rather  than  being 


inherently  specific  to  a  single  flow  condition. 


We  also  emphasize  the  effects  of  the  added  mass  associated  with  submerged  coating 


response.  In  addition,  we  have  extended  our  previous  analyses  in  order  to  consider 


three-dimensional  flow  distributions  actively  forcing  three  dimensional  coating  surfaces 


and  we  have  initiated  study  of  direct,  interactive  dynamic  coupling  between  the 


responding  surface  and  a  surface-altered  flow  field.  All  of  the  calculations  described 


here  were  processed  on  CRAY  I  computers. 


V.v.v. 


Initial  simulations  of  material  response  arc  non- interactive  between  fluid  and  solid 

regimes.  Instead,  the  surface  excitations  arc  independently  forced  by  modeled  random 

space/time  turbulent  pressure  distributions  to  obtain  information  on  material  response 

characteristics.  These  surface  excitation  simulations,  mimic  a  turbulent  pressure  field  on 

a  non-compliant  flat  plate.  These  are  generated  by  application  of  the  Ash-Khorrami 
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Monte- Carlo  turbulent  pressure  field  model.  '  These  computations,  albeit 

non- interactive,  are  useful  for  screening  relative  response  influence  of  materials, 

material  properties,  composites,  and  internally  structured  coating  candidates.  We  thereby 

reduce  the  number  of  candidate  coatings  and  flow  conditions  that  need  to  be  considered  in 

subsequent  fully  interactive  fluid-solid  coupling  simulations.  The  utility  of 

non-interactive  procedures  in  this  screening  process  is  particularly  evident  in  analyses  of 

the  isolated  influences  of  overburden  in  ocean  depth  studies,  virtual  mass-surface 

displacement  modulation,  and  three-dimensional  influences  relative  to  the  theoretically 

more  familiar  two-dimensional,  plane  strain  influences. 

For  direct  fluid/solid  interactions,  Navier-Stokes  flow  simulations  are  used.  We 

6  7 

apply  a  pseudo-spectral  method  based  on  Orszag  and  Kells  and  Orszag  and  Patera,  later 

g 

extended  to  include  small  displacement  rate  influences  of  linearized  surface  response. 
Flow  instabilities  of  the  Tollmein-Schlichting  form  are  impressed  on  otherwise  stable 
viscous  profiles.  The  influence  of  detailed  solid  coating  dynamic  compliance  is  traced  in 
the  subsequent  growth  or  reduction  of  the  viscous  instabilities. 

II.  Simulations  &  Conditions 

Finite  Element.  Procedures 

For  the  compliant  coating  materials  response  simulations  we  apply  a  discretized 
procedure  to  simulate  the  response  of  composite  materials  and  internally  structured 
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coatings.  To  follow  high  frequency  and  relatively  large  deformation  motions,  we  adapt 

both  implicit  and  explicit  two-dimensional  and  three-dimensional  versions  of  the  Hallquist 
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finite  element  codes.  '  In  some  of  our  computations,  explicit-in-time  finite  element 
computations  are  applied  since  the  Courant-stable  time  interval  for  elemental  response  is 
of  the  same  order  (or  larger)  than  the  virtual,  implicit  time  interval  sampling  rate 
required  for  the  response  frequency  range  of  interest.  Having  established  the  explicit 
fluid/solid  time  data  base  associated  with  the  dynamic  response  of  representative 
coatings,  we  are  able  to  develop  interpolation  and  time-sequencing  (sub  cycling)  rules  and 
algorithms  for  automatic  fluid/solid  coupling  procedures  in  subsequent,  efficient  implicit 
code  applications. 

Discussion  of  the  finite  element  procedure  is  based  on  the  implicit,  iterative 
formulation,  for  generality.  The  unsteady  dynamical  problem  is  represented  as  an 
incremental  (A)  perturbation  to  the  equilibrium  state.  A  lumped  or  consistent  mass 
matrix  relationship  describes  the  instantaneous  perturbation  at  an  incremental  time 
epoch,  N  +  1,  based  on  a  solution  (or  initial  conditions)  at  a  previous  time  epoch,  N. 


M?N+1)+  D(d(N)  +  oAd*(N))d(N+1)+  K^X(N))  A^, 


=  p(x(N))(N+1)  _  F  (X(N))  , 


where  M  is  the  lumped  mass  matrix,  d  is  the  nodal  displacement  vector  relative  to 

coordinate  vector,  d  =  (X^+^  -  X^),  is  the  material  property  dependent,  tangent 

stiffness  matrix,  Ad*v  is  the  iteratively  estimated  displacement  increment  for  this  load 

step  while  Ad^  is  the  prescribed  displacement  increment  for  the  implicit  time  increment, 
N+l  N 

At  =  t  -  t  .  P  is  the  external  loads  vector  based  on  the  applied  loading  at  epoch  N+l 
but  on  the  geometry  at  epoch  N.  The  stress  divergence  vector.  F,  is  based  on  the 


displaced  state  and  stresses  at  load  step  N.  Time  rates  of  change  and  accelerations  of  the 
displacement  vector  are  symbolized  by  the  single  and  double  dot  superscripts,  respectively. 

The  Rayleigh  damping  matrix  is  iteratively  evaluated  as  linear  combination  of  mass 
and  stiffness  matrices. 

D  =  aM  +  BKt  ,  (2) 

where  a  and  g  are  user-prescribed  constant  scale  factors  with  a  typically  less  than  or 
equal  to  1/2  and  g  typically  equal  to  1/4.  Time  evolution  of  equation  (1)  involves  use  of  an 
unconditionally  stable,  one  step,  Newmark-g  time  integration  scheme  in  conjunction  with 
a  choice  between  one  of  three  (2  quasi-Newton  an  one  modified  Newton- Raphson) 
iterative  matrix  reduction  procedures,  based  on  the  Green-Nagdhi  stress  rate 
formulation.  The  method  is  stable  for  a  broad  range  of  nonlinear  probelms  involving  finite 
strains  and  (arbitrarily)  large  rotations.  A  line  search  is  used  with  the  iterative  step  that 
is  combined  with  automatic  matrix  stiffness  reformation  to  avoid  nonconvergence. 

Initial  and  boundary  condition  prescriptions  are  flexible.  Virtual  mass,  inertial 
influences,  localized  nodal  loading,  exterior  pressure  distributions  and  time  histories, 
surface  and  internal  shear  boundary  conditions,  internal  temperature  sources, 
distributions,  and  conduction  and  thermal  material  property  changes  are  all  provided 
together  with  a  self-consistent  Jaumann  stress  rate  formulation  and  a  variety  of 
constitutive  relations,  for  generality.  In  addition,  slide  line  logic  (or  slide  surface  logic  in 
the  three  dimensional  studies)  permits  treatment  of  discontinuous  internal  material  layers 
and  imbedded  material  discontinuities,  in  these  simulations. 

Concentrated  displacement  inertial  nodal  mass  loading  permits  automatic  evaluation 
of  added  mass  effects  but  with  some  penalty  in  computing  efficiency,  since  the 


consistent  mass  matrix  formulation  requires  incremental  re-evaluation  at  every  stage  of 
elemental  displacement.  Since  we  are  also  interested  in  studying  the  effects  of 
submergence  infuences  (overburden  and  temperature)  on  material  response,  we  choose  a 
simpler,  less  time-consuming,  but  equally  valid  approach.  We  simply  apply  a  "thin*’  layer 
of  water  over  the  coating  surface  to  generate  shallow,  depth  added  mass  inertial  influence 
on  surface  deflection  response.  "Thin"  here  refers  to  the  relative  fractional  depth  of  this 
water  overlayer  with  respect  to  the  underlying  depth  of  coating  material  (typically  1/5  or 
less).  This  provides  adequate  "coverage"  to  model  coating  surface  deflection  added  mass 
inertial  influences  without  incremental  re-evalution  of  the  mass  matrix.  Of  course  in  our 
depth  of  submergence  studies,  discussed  later,  the  layer  of  ocean  water  covering  the 
coating  was  deepened  appropriately  for  imposition  of  consistent,  static  overburden  as  well 
as  added  mass  influences. 

Material  Models  and  Flow  Condition 

Of  the  six  general  classes  of  material  models  available  in  the  finite  element 

codes,  our  analyses  here  are  all  based  on  the  use  of  only  two:  elastic  and  linear 

visco-elastic.  Fig.  1  illustrates  the  general  flow  and  material  model  geometry.  While  we 
12 

have  examined  single  layer  and  multi-layered  combination  compliant  coatings 
consisting  of  natural  rubber,  foams,  RTV  silicone  rubber,  polyurethane  and  other 
materials,  we  are  presently  concentrating  on  promising  combinations  of  basically  two 
materials.  The  softest,  most  yielding  of  these  is  polyvinylchloride  (PVC).  This  is  a 
time-cured,  gelatinous,  visco-elastic  material  whose  shear  modulus  may  vary  by  one  or 
more  orders  of  magnitude  with  small  changes  (a  few  percent)  in  the  ratio  of  the  polymer 
to  its  binder  and  filler.  The  other,  stiffer  material  is  the  more  nearly  elastic  Neoprene, 
which  is  conceptually  used  as  an  overlayer  and  for  supporting  internal  structures  in  the 
compliant  coating  candidates  discussed  here. 


Figure  1  represents  a  modeled  planar  panel  of  coating  material  flush  mounted  in  a 
semi-infinite  flat  plate  (symbolized  by  the  dashed  lines  extending  indefinitely 
streamwise).  The  cartesian  coordinate  orientation  for  two-dimensional  (x,z)  and 
three-dimensional  (x,y,z)  computations  are  sketched  at  the  top  of  the  figure  together  with 
the  mean  potential  flow,  U  ,  orientation.  Streamwise  coating  lengths  ranging  from  12  cm 
to  60  cm  and  spanwise  widths  from  5  cm  to  12  cm  (for  three-dimensional  computations) 
have  been  considered.  Total  coating  depths  of  from  2  to  8  cm  have  been  studied.  Panel 
locations  for  the  application  of  uncoupled  (surface  excitation)  or  fully  fluid/solid 
interactively  coupled  pressures  establish  the  finite  element  streamwise  and  spanwise 
nodal  spatial  distribution.  From  10  to  64  streamwise  elements  and  from  3  to  5  spanwise 
elements  have  been  used  while  the  in-<iepth  distribution  of  elemental  spacing  varies  from 
3  to  12  elemental  segments  in  our  computations.  The  two-dimensional  plane-strain 
computations  simulate  a  vertical  cut  in  a  semi-infinite  span  of  coating  materials  while 
the  three-dimensional  computations  incorporate  finite  span  and  spanwise  edge  influences. 
Streamwise  and  spanwise  edge  boundary  conditions  tested  included  periodic,  fixed 
(adherent),  or  vertically  free  (end  roller)  designations.  The  base  of  the  coating  panel  was 
modeled  with  a  rigid  constraint,  when  it  was  found  that  the  modeled  results  were  identical 
using  either  this  simple  designation  or  an  actual  nearly  rigid  metal  surface  (steel) 
elemental  substrate  layer. 

We  have  computed  results  for  mean  flows  of:  1.0,  1.52,  3.0,  10.0,  and  30.0  m/s  at 
depths  from  "ocean  surface"  (submerged  to  a  depth  of  about  1/5  to  1  coating  depth 
thickness)  down  to  200  m  with  static  depth  overburden  and  temperature  variations 
included  as  initial/boundary  conditions. 

Two  classes  of  compliant  coatings,  each  combining  both  materials,  Neoprene  and 
PVC,  are  studied  here.  One  consists  of  a  thin,  homogeneous  overlayer  of  the  stiffer 
Neoprene  (shown  darkened  in  Fig.l)  over  the  softer  PVC  underlayer.  The  Neoprene 
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overlayer  thickness  in  our  calculations  varied  from  1/5  to  1/10  the  thickness  of  the 
underlaying  PVC.  The  second  class  of  compliant  coatings  also  uses  the  thin  Neoprene 
overlayer  but  includes  internal  structure  consisting  of  spanwise  rectangular  slots  with  a 
slot  depth  (or  height)  about  1/3  of  the  PVC  layer  thickness  and  breadth  or  width  varying 
over  range  varying  from  1/2  to  4  times  the  slot  depth.  Neoprene  sidewalls,  conceptually 
separated  and  supported  the  unbedded  slots.  Computations  included  considerations  of 
slots  filled  with  water,  with  PVC  gel,  or  empty  slots  (voids). 

The  viscoelastic  shear  modulus  of  the  Neoprene  and  PVC  are  modeled  by  a  standard 
linear  time  (t)  relaxation  procedure  for  linear  solids^ 1  represented  by, 

G(t)  =  G  +(G  -G  )e"Bt  (3) 

00  VJ  00 

in  which  G(t)  is  the  composite  shear  modulus,  while  and  Gq  are  long  and  short  interval 

shear  moduli,  respectively,  and  B  is  the  relaxation  frequency  parameter. 

Figure  2  illustrates  the  experimentally  based  (and  analytically  fit)  dependence  of  the 

shear  modulus  on  frequency  for  the  relatively  stiff  Neoprene  at  a  temperature  of  at  10°C, 

and  for  the  much  more  compliant  PVC  at  10°C  and  softer  yet  at  20°C.  The  dependence 

13 

shown  closely  correlates  the  experimental  values  reported  by  Hunston  and 
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Gad-El-Hak  .  The  density  of  the  PVC  varied,  in  our  model,  from  500  to  1025  Kg/m  , 

depending  on  the  polymeric  percentage  in  the  composition  mix.  The  Neoprene  mass 

3 

density  was  assumed  constant  at  1250  Kg/m  .  Both  Neoprene  and  PVC  are  essentially 
incompressible  with  a  poisson  ratio,  v=l/2.  The  bulk  moduli  of  the  PVC  and  Neoprene 
were  set  at  620  and  5270  MPa  respectively,  in  the  model. 

Screening  Non-interactive  "Monte- Carlo”  Pressure  Exicitation 

Materials,  composites  and  internally  structured  coatings  need  to  be  screened  and 
reduced  in  number  so  that  more  exacting  coupled  computer  simulations  and  experiments 


can  be  limited  to  a  practical  number.  Our  screening  procedure  makes  use  of  an  existing 
Monte-Carlo  turbulent  pressure  excitation  computer  model  that  can  be  "patched"  for 
boundary  conditions  to  detailed  finite  element  representations  of  the  coating  response. 

4  5 

Ash  and  Ash  and  Khorrami  have  developed  a  Monte-Carlo  procedure  for  modeling 
turbulence  pressure  excitations  on  a  flat  plate.  Statistical  turbulent  functions  are 
simulated  that  are  consistent  and  convergent  to  within  arbitrarily  close  representation  of 
experimental  measurements . 1 

To  initiate  a  calculation  of  the  turbulent  pressure  we  consider  the  flow  geometry 
over  a  flat  plate  with  zero  pressure  gradient,  immersed  immediately  below  the  surface  in 
sea  water.  We  examine  conditions  over  a  range  of  specified  steady  state  free  stream 
velocities.  The  turbulent  field  development  is  completed  following  transition  over  a 
specified  model  transitional  length.  Also  specified  is  the  minimum  turbulent  structure 
wavelength  of  interest.  The  Ash  code  generates  only  a  random  space  and  time  varying 
pressure.  Convection  and  convection-dominated  burst-streak  growth  and  upthrust 
vorticity  boundary  layer  features  are  separately  imposed  as  model  parameters. 

Figure  3  shows  the  Ash  model  results  on  a  Fourier  transformed  pressure  frequency 
spectra,  for  a  modeled  10  m/s  ocean  water  flow  at  a  spatial  position  25%  from  the  leading 
edge  of  flat  plate  which  is  0.32  m  long.  The  calculation  used  a  specified  transition  length, 
upstream  of  the  model  compliant  panel  leading  edge,  of  0.50  m. 

Time  and  space  scales  to  be  resolved  with  the  pressure  model  calculations  are 
dictated  by  the  desired  structural  frequency  response  range  we  wish  to  examine.  This 
range  is  controlled  by  finite  elemental  nodel  spacing,  material  properties  of  the  coating, 
the  presence  or  absence  of  voids,  and  the  overall  panel  dimensions. 

Once  the  pressure  patch  statistics,  particularly  the  power  spectra  and  decay 
functions,  are  checked  and  found  to  be  reasonable  representations  of  experimental 


pressure  distribution  histories  these  statistics  are  saved  on  disk.  The  same  excitation 
model  is  used  to  stimulate  the  response  of  many  selected  coating  panel  candidates  using 
identical  space  and  time  scales. 

Interactive  Fluid/Coating  Coupling  With  the  Full  Navier-Stokes  Equations 

In  this  phase  of  our  analysis,  the  solid  material  response  calculations  are  simulated 
using  the  Hallquist  2- and  3-dimensional  finite  element  computer  code  methods.9,1^'11 
Initially  we  are  concentrating  our  efforts  on  examining  material  coatings  that  may  delay 
transition  to  turbulence  by  reducing  the  effects  of  laminar  instabilities  at  sub-transitional 
Reynolds  numbers. 

Our  current  studies  impose  single  wave  Tollmein-Schlichting  instabilities  on  a 
Blasius  flow  profile,  extending  from  a  level  coating  surface  to  a  mid-plane  region  that  is 
outside  of  the  significant  viscous  boundary  layer  influences. 

Our  primary  consieration  here  is  to  develop  and  apply  a  method  for  efficiently  and 
interactively  coupling  a  pseudo-spectral  Navier-Stokes  simulation  of  the  unsteady, 
transitional  boundary  layer  flow  past  a  compliant  surface  represented  by  finite-element 

g 

description.  Our  fluid  dynamic  procedure  uses  a  modification  of  the  FLOGUN  code, 

6  7 

which  is  based  on  the  channel  flow  procedures  by  Orszag  and  his  colleagues.  Surface 
displacement  rate  is  imposed  as  a  boundary  condition  for  the  velocity  of  the  fluid  at  the 
wall.  These  displacement  rates  are  passed  to  the  fluid  case  incrementally  from  the  finite 
element  code.  The  fluid  code  computation,  in  turn,  is  used  to  determine  the  cyclic 
variations  in  pressure  distribution  that  drive  the  wall.  The  interactive  coupling  proceeds 
indefinitely  using  either  direct  or  implicit  time-step  increments.  The  equations 
describing  the  flow  field  are  the  incompressible  Navier-Stokes  equations. 


and  the  incompressibility  condition, 


V»v(x,t)  =  0 


where  v  is  the  velocity,  and  x  is  the  position  in  two-dimensional  space.  We  have  also 
introduced  p,  the  pressure,  (divided  by  the  density),  and  v,  the  kinematic  viscosity.  The 
domain  of  integration  is  given  by 


0  <  x  <  X  , 


0  <  z  <  oo  , 


with  periodic  boundary  conditions  in  the  x-direction  (the  horizontal  direction)  and 
free-stream  inviscid,  potential  flow  boundary  conditions  at  z  =  <».  We  impose  the 
free-stream  boundary  conditions  at  a  normal  distance  from  the  wall  approximately  five 
times  the  local  boundary  layer  thickness.  At  the  wall,  the  tangential  velocity  is  zero, 
while  the  normal  velocity  is  equal  to  the  displacement  rates  given  by  the  finite  element 
results  for  the  wall  motion.  For  simplicity,  this  boundary  condition  is  imposed  at  z  =  0. 
The  velocity  is  considered  a  smoothly  varying  function  and  is  represented  by  a  truncated 
series  expansion  in  spectral  space. 


v(x. t)  =  £  I  u(m,p,t)  exp[2Trimx/X]  T  (z)  .  (7) 

|  m  |  <M  p=0  p 

Here  Tp(z)  -  cos  (p  cos_1z),  the  Chebyshev  polynomial  of  degree  p.  Using  Chebyshev 
polynomials  in  the  direction  normal  to  the  wall  results  in  good  resolution  in  the  boundary 
layer,  with  about  one  third  of  the  collocation  points  being  clustered  there.  Efficiency  is 
gained  by  evaluating  spatial  derivatives  in  spectral  space  then  transforming  back  to 
physical  space  for  evaluation  of  the  non-linear  advective  products.  Since  the  Chebyshev 
expansion  is  equivalent  to  Fourier  expansion  on  a  distorted  space,  fast  Fourier  transforms 
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can  be  applied  either  to  the  normal- to-  the-wall  transformation  between  Chebyshev  space 
and  physical  space  or  the  parallel-to-  the-wall  transformation  between  Fourier  space  and 
physical  space. 

Time  evolution  is  accomplished  in  three  steps  with,  for  example,  the  vorticity 
evaluated  in  spectral  space,  then  transformed  to  physical  space.  The  vxu  contribution  is 
evaluated  in  physical  space.  This  step  is  performed  explicitly  in  time,  with  an 
Adams-Bashforth  method. 

Development  of  specialized  velocity/pressure  boundary  conditions  and  algorithums 
for  the  coupled  fluid/solid  interface  flow  are  described  separately  in  a  current 
publication,  ^  so  that  we  will  not  dwell  on  them  here.  Results  of  these  coupled  studies 
indicate  that  a  substantial  amount  of  the  disturbance  wave  energy  is  transferred, 
in-phase,  to  the  solid  where  it  is  absorbed  in  resonant,  displacement-convection  process, 
in  phase  with  the  advected  flow  disturbance.  This  flow-solid  boundary  interaction  event 
is  in  contrast  to  coincident  fluid/solid  dispersion  processes  which  are  justifiably  receiving 
considerable  attention.8,13'14  The  coupled  simulations,  of  course,  involve  both 
convective  and  dispersive  (as  well  as  static)  influences  generally.  It  remains  an  important 
part  of  our  task  to  unravel  and  identify  these  separate  effects  in  our  present  and 
continuing  analysis. 

III.  Discussion  of  Results 

This  discussion  emphasizes  new  and  different  results  and  suggestions  about 
compliant  coating  candidates  in  comparison  to  our  previous  summary.  In  particular,  we 
address  here:  the  special  influences  of  fully  three  dimensional  geometries  vs  two 
dimensional  plane  strain  results,  the  significant  influences  of  added  mass  and  depth 
overburden  in  comparison  to  our  previous  in-vacuo  results,  multicomponent  layered 


coatings  including  the  ordered  sequence  of  materials  and  layer  thickness,  and  internally 
structured  coatings  including  considerations  on  aspect  ratio  of  the  cavities  and  the 
influence  of  cavity  fillers. 


Three  Dimensional  vs.  Two  Dimensional  Results 

The  spanwise  edges  of  the  compliant  surface  panel  may  be  considered  partially  or 
fully  constrained  within  a  flat  plate  cavity  in  analogy  with  a  physical  flow  experiment. 
The  constraint,  or  lack  of  it,  alters  the  character  of  the  disturbance  waves  that  are 
reflected  from  these  boundaries,  and  significantly  alters  the  response  time  histories  of  the 
coating  candidates  from  the  values  developed  using  a  strictly  two-dimensional,  plane 
strain  analysis.  In  addition,  the  boundary  layer  flow  itself  possesses  substantial  spanwise 
inhomogenie  ties  in  flow  structure. 

4  5 

Using  the  non-interactive  Ash-  Khorrami  ’  Monte-Carlo  turbulent  pressure  model, 
a  one-dimensional  streamwise  space/time  varying  pressure  "strip"  is  generated  for 
exciting  the  two  dimensional  (finite  length  and  finite  depth,  infinite  span)  plane  strain 


results.  For  our  three  dimensional  panel,  parallel  "strips"  are  applied  along  adjacent  rows 
of  elements.  The  adjacent  strips  are  considered  to  be  statistically  independent,  that  is, 
spanwise  uncorrelated.  This  is  physically  analogous  to  assuming  that  each  row  of 
elements  is  populated  with  surface  microphones  along  its  length  (the  pressure  sensitive 
surface  "nodes"  of  the  individual  finite  elements),  yet  laterally  separated  at  a  modest  but 
finite  distance  from  each  other  that  is  greater  than  the  spanwise  correlation  length  scale 
of  the  turbulent  pressure  field  represented. 

To  illustrate  the  general  influences  of  the  three-dimensional  panel  results  we  first 
present  two-dimensional  surface  deformation  response  spectra  in  Fig  4  for  a  10  m/s  ocean 
water  flow  in  zero  mean  pressure  gradient  over  a  3  cm  thick  panel  of  polyvinylchloride 
(PVC). 
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Three-dimensional  results  for  this  flow  condition  and  coating,  using  spanwise 
unconstrained  (free  surface  edge)  panel  edge  conditions,  is  shown  in  Fig.  5.  Volumetric 
propagation  and  reflection  of  the  disturbance  waves  from  finite  span  extent  together  with 
three-dimensional  material  deformations  are  seen  to  modulate  and  interactively  reduce 
the  deformation  amplitude  and  response  frequency  content  at  the  higher  frequencies 
(much  as  a  low  pass  filter  applied  to  a  partially  damped  inertial  system). 

In  practice,  the  two-dimensional  plane  strain  simulations  are  used  for  the  primary 
part  of  the  material  model  and  structural  screening.  The  three-dimensional  results 
continue  to  be  reserved  for  estimates  on  three  dimensional  geometrical  and  flow  structure 
influences  on  a  limited  number  of  the  more  promising  candidates  because  of  their  much 
greater  cost  and  analytical  complexity. 

Overburden  and  Added  Mass 

We  tested  and  isolated  the  influence  of  added  mass  (water  inertia)  on  our  surface 
response  computations  and  evaluated  the  results  in  comparison  to  in-vacuo  response 
results  for  a  10  m/s  ocean  water  flow.  We  compare  in-vacuo  displacement  spectra.  Fig. 
6,  with  added  mass  results  in  Fig.  7.  These  plane  strain,  2D  tests  were  computed  for  the 
midplane  position  on  a  test  panel  of  PVC  assumed  to  be  1  cm  thick.  The  added  mass 
effects  reduce  the  response  to  negligible  levels  in  the  intermediate  range,  60  to  100  Hz, 
and  somewhat  modulate  the  results  at  even  lower  frequencies.  At  higher  frequencies  yet, 
the  influence  of  virtual  mass  is  not  influential  and  at  about  one  kilohertz  is  undetectable, 
as  linear  resonance  analysis  suggests.  However,  since  our  attention  here  is  specifically  on 
the  frequency  band  from  1  Hz  to  about  200  Hz  surface  response,  we  include  added  mass 
influences  in  all  of  our  current  and  future  simulations  on  the  surface  coating  response. 

Another  important  influence  that  must  be  considered  in  analysis  of  compliant 
coatings  for  submersible  applications  is  the  effect  of  ocean  depth.  Unlike  the  added  mass 


inertial  correction,  submergence  to  great  depth  with  a  compliant  coating  over  a  rigid 
substrate  effectively  pre-stresses  the  material  to  a  level  necessary  to  equilibrate  the 
hydrostatic  overburden  of  the  water.  In  effect,  the  same  hydrostatic  overburden  that 
applies  a  crushing  force  to  submersible  hulls,  stiffens  the  softer  compliant  coatings  with 
pre-stress  loadings  that  approach  or  exceed  their  elastic  moduli.  The  test  calculations 
here  were  initiated  at  an  internal  material  stress  state  equivalent  to  the  overburden 
(depth)  condition  of  interest.  Amplitude,  wavelength  and  frequency  shifts  were  computed 
as  a  function  of  overburden.  The  influence  of  temperature  (cooling)  at  ocean  depth  on 
material  properties  was  also  tested  but  found  to  be  negligible  with  respect  to  the 
overburden  influences. 

The  overburden  stiffening  with  depth  was  found  to  result  in  a  monotone  decreasing 
value  of  the  root  mean  square  amplitude  of  coating  normal  deflection  ratio  (A/A^)^  for 
both  Neoprene  and  softer  PVC  coatings.  As  anticipated,  the  effect  of  depth  is  greater  for 
the  softer,  lower  elastic  modulus  material,  PVC,  than  for  the  stiff er  coating.  Neoprene,  as 
seen  in  Fig.  8.  The  influence  of  depth  on  the  band  width  of  a  significant  portion  of  the 
frequency  response  of  PVC  coating  is  shown  in  Fig.  9.  The  bandwidth  here  is  that 
frequency  range  of  the  deformation  spectrum  in  which  90%  of  the  significant  deformation 
occurs.  The  corresponding  influences  on  the  much  stiffer  material.  Neoprene,  were  so 
insignificant  that  they  were  not  plotted  here. 

Layered  Homogeneous  Coatings 

The  thickness  of  the  compliant  coating  is  a  significant  factor  in  response  amplitude. 

whether  it  appears  as  a  single  layer  or  part  of  a  multiply  layered  coating  over  a  stiff 

substrate.  Increasing  thickness,  however,  is  not,  by  itself,  a  positive  influence  on 

desirable  compliant  coating  response,  since  a  too-thick  depth  of  coating  can  be  a  factor 

13,14 


leading  to  the  onset  of  static  divergence. 


The  resulting  larger  amplitude  surface 


deformations  generally  occur  at  very  low  frequency,  which,  as  in  the  static  divergence 
case,  may  offset  drag  reduction,  and  may,  in  fact,  increase  the  drag. 

Our  results,  however,  show  that  use  of  ordered  multilayer  coatings  rnay  produce  the 
sought  for  results  on  controlled  small  amplitude  deformation  at  substantial  frequencies. 
This  is  particularly  true  with  a  thick  (1  cm  to  2  cm)  layer  of  soft  material  (PVC)  over 
which  is  placed  a  much  thinner  (2  mm)  coating  of  a  stiffer  material  (Neoprene).  The 
results  at  a  mean  flow  speed  of  about  10  m/s  appear  encouraging  using  such  composite 
multilayer  coatings  in  our  simulations. 

Internally  Structured  Coatings 

Our  analysis  on  internally  structured  coatings,  to  date,  has  been  restricted  to 
coatings  with  internal  reinforcing  stiffeners  and  channeled  openings  that  are  oriented 
transverse  (spanwise)  to  the  mean  flow  direction  in  the  coating.  We  have  examined  both 
liquid  filled  and  empty  (void)  channels  using  a  viscous  liquid  (PVC)  and  an  incompressible 
liquid  (water)  as  a  filler  in  comparison  to  open,  coating-imbedded,  channels  (voids).  The 
voids  yield  the  most  promising  results  in  terms  of  rms  and  peak  response  amplitudes  to 
surface  pressure  excitation. 

Shallow  but  wide  channels,  particularly  for  voided  (no  filler)  cases,  with  channel 
width  to  height  of  3  or  more  also  show  the  most  promise  for  coating  response  amplitude 
enhancement  as  seen  in  Fig.  10.  Our  simulations  incorporate  the  thin  outer  stiffened 
layer  of  firmer  material  as  an  upper  layer  to  the  stub  supported  channels  within  a  thick 
layer  of  the  low  modulus  PVC  material.  Excitation  induced  frequency  response  of  the 
structured  coatings  is  also  much  more  satisfactorily  reproduced  over  a  broad  range  of 
excitation  frequencies  when  voids  rather  than  liquid  filler  are  used. 
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IV.  Summary 

Some  preliminary  indications  apparent  from  our  numerical  studies  at  this  time,  arc 
summarized. 

•  Added  mass  and  depth  of  submergence  (overburden)  considerations  must  be  included 
in  analysis  of  compliant  coating  performance,  particularly  for  the  softest,  lowest 
moduli  materials,  which  dynamically  respond  as  if  they  were  much  stiffer  materials 
when  those  influences  are  included  in  our  computations. 

•  Composite  multilayered  coatings  show  promise  for  drag  reduction,  particularly  when 
thickened  soft  inner  coating  layers  are  sandwiched  between  a  stiffer  thin  outer 
coating  layer  and  a  rigid  substrate. 

•  Shallow  flat  channels  imbedded  in  the  material  and  oriented  transversely  (spanwise) 
to  the  mean  flow  direction  may  favorably  modulate  amplitude,  frequency  and 
wavelength  of  the  response,  particularly  if  the  channel  is  left  open,  free  of  filler 
material. 

•  Three  dimensional  influences  of  both  the  flow  patterns  and  the  compliant  surface 
geometry  should  be  evaluated  and  corrected  for,  prior  to  recommending  promising 
coating  candidates  for  further  experimental  verification. 
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FIGURE  CAPTIONS 


Fig.  1.  Geometry  of  compliant  surfaces,  flow  and  coordinate  orientation  for  finite 
element  transient  response  simulations. 


Fig.  2.  Shear  modulus  of  two  critical  materials  as  a  function  of  frequency  for 


compliant  coating  response  models. 


Fig.  3.  Driving  pressure  excitation  spectral  distribution  from  the  Ash  Monte-Carlo 
model  (Refs.  4  and  5)  for  a  10  m/s  ocean  water  flow  at  a  position  25%  from  the  leading 
edge  of  a  0.32  m  flat  plate  compliant  surface. 


Fig.  4.  Fourier  transformed  displacement  vs  frequency  distribution  for  the  response  at 
the  mid-panel  position  of  a  3  cm  thick  PVC  panel  in  two-dimensional  flow  at  10  m/s. 


Fig.  5.  Fourier  transformed  displacement  vs  frequency  distribution  for  the  response  at 
the  mid-panel  position  of  a  3  cm  thick  PVC  panel  in  three-dimensional  flow  at  10  m/s. 
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Fig.  6.  Fourier  transformed  displacement  vs  frequency  distribution  for  the  response  at 
the  mid-panel  position  of  a  1  cm  thick  PVC  panel  in  a  two-dimensional  flow  at  10  m/s. 


Fig.  7.  Fourier  transformed  displacement  vs  frequency  distribution  for  the  response  at 
the  mid-panel  position  of  a  1  cm  thick  PVC  panel  including  added  mass  in  a 
two-  dimensional  flow  at  10  m/s. 
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